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Abstract: The stand growth and yield dynamic models for Larch in Jilin Province were developed based on the forest growth 
theories with the forest continuous inventory data. The results indicated that the developed models had high precision, and they 
could be used for the updating data of inventory of planning and designing and optimal decision of forest management. 
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Introduction 

Forest growth and yield models have been studied 
widely and intensively. For example, basal area models 
were developed for slash pine plantations by Pienaar et al. 
(1981), and for loblolly pine plantations by Sulliyan et al. 
(1972) and Burkhart et al. (1984), the survivor tree models 
were constructed for loblolly pine by Lemin et al. (1983), 
for radiata pine by Pienaar et al. (1973), for slash pine by 
Bailey et al. (1985), and both basal area and survivor tree 
stem models for loblolly pine were developed by Cao et al. 
(1982) 

In China, however, the research on growth models of 
forest has been limited due to the fact that data sets used 
in analysis were from temporary plots and that there was 
no enough permanent plot data,. With the development of 
continuous forest inventory (CFI), tens of thousands of 
permanent plots have been measured and re-measured, 
which provides very high quality information on many fac¬ 
tors. These valuable data of the plots warrant further de¬ 
velopment and application forest growth and yield models. 
The objective of this paper is to use CFI data for develop¬ 
ing dynamic growth and yield models for larch in Jilin 
Province in China. The most recent inventory in Jilin indi¬ 
cated that larch comprises a high proportion of the forests. 

Data collection 

The data used in this study came from the 
re-measurements of Larch plantations in the CFI in 1984 
and 1989 in eastern Jilin. Plot is 0.06 hm 2 in size. The di¬ 
ameter at breath height (DBH) of every tree in the plot and 
all stand factors are measured except the dominant height. 
For this research, only plots on which the trees had been 


measured on both occasions, the basal area was more 
than 65% of the whole stand and the crown closure was 
above 0.5, were used. The number of available plots was 
112, and 82 of them were used for setting up the models 
and 30 for testing the models. 

Development and selection of models 

Mathematical expressions for static growth and yield 
models are as follows: 

For a particular site quality; Y=f(S, A) 

For variable site quality: V=/(S, A, SQ) 
where S is stand density; A is stand age; SQ is site quality; 
and Y is stand factors. 

The feature of this kind of models is that the predictors 
are independent of the stand state at the beginning of the 
prediction. 

Dynamic models are those whose predictors relate to 
the state of the stand at the beginning of the prediction. 
They can be expressed as: 

For a particular site quality: Y z =f{S, Ai, A 2 , yi) 

For variable site quality: Y 2 =f(S, A u A 2 , SQ, yj) 
where Y 2 is prediction of the stand factors at A 2 of stand 
age; Yi is measurement of the stand factors at Ai of stand 
age, and A Z ^A : . It is obvious that Y 2 should be equal to X 
when A 2 equals Ai. 

After some experimentation with the following models 
(equations 1 to 3), we found that these models were not 
very suitable for larch. Thus some more flexible models 
were needed. 


ln(X 2 ) = 


—xln(X, )+ax 
A, 



+ bx 



(D 


where a, b are parameters. 
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ln( Y 2 ) = a 0 + a, x SI 


-x (a 2 + a 3 x SI + ln(y,)) 


where ao, ai, a 2 , and a 3 are parameters. 

Researches by Pienaar et at. (1973), Weng (1989), and 
Li (1990) indicated that the parameter k could be described 
as a function of stand density and parameter a could be 
described as a function of site quality in Richards equation: 


Y = ax(l-e ( - KxA) ) c (4) 

where K and c are parameters. 

That is, equation (4) could be described as: 

Y = f(SQ)x(l- e , “ /(4)x ' 4) ) 1 (5) 

It can be assumed that the growth of stand factors is A Y 
when the length of the growth period is A A. If A A is small 
enough, then 


AY = f{SQ)x{l-e 


l-f(S)xAA) V 


Equation (6) is a function of growth. But it is a monotonic 
function, thus it is not reasonable theoretically. Especially it 
may produce large errors for fast growing species. There¬ 
fore some modification is needed for equation (6). For 
even-aged mono-species stands, stand growth depends on 
not only the length of period, but also stand age. Thus the 
term exp(d*A ; ) is used as a modification factor. Equation (7) 
can then be obtained: 


AY = f(SQ)x(\-e 


A J — A 2 + 


(d - c x e 1 1 <SI> ) 


Estimation of parameters and test and use of the 
models 

Estimation of model parameters 

Evaluation of site quality 

Although dominant height was not measured, a site 
class table can be constructed to evaluate site quality with 
average stand height. The average stand height is defined 
as average height of all trees on the plot. To integrate site 
quality into the models, the average stand height at base 
age was used as an index of site quality. This index is dif¬ 
ferent from site index, which is referred as site class index. 
The equation of site class index is: 


(-.02035 xA„ I , 
_<-.02035 xA), 


where SCI is site class index; Ao is base age (30 years); A 
is stand age; H is average stand height. 

Determinations of the site quality function and stand den¬ 
sity function 

For simplicity, the site quality function is defined as: 


(-f(S)xAA) 


f (SQ ) = a'xS/ 


and the stand density function is described as 


f(S) = kyY { or f(S)=ky\ og$;) 


where d is parameter. 

Because A'i=Y 2 -Y, and A A=A^A t , they can be sub¬ 
stituted into equation (7), and the final dynamic growth 
model is: 

Y 2 = Y i + f (SQ)x(\ - e < ~ J(SMA -' A,)) f xe ,dyA,> ' (8) 

Equation (8) is not a monotonic function and has only 
one maximum value. The beginning age of growth period 
at the maximum point is as follows: 


where a’, b and k’ are parameters. 

Results 

Equation (8) can be solved by least square. Then the op¬ 
timal parameters are estimated with Maquardt logarithm 
using the least square parameter estimates. The goodness 
of fit criteria is residual sums of square and the coefficient 
of determination. 

The results indicate that equation (8) is much better than 
equations (1), (2), (3), and (6). The parameters and the 
goodness of fit criteria are list in Table 1. 


Table 1. Parameters and Goodness of Fit of Models 

Stand Factors _Density Function_Parameters_ Goodness of fit 




a' 

b 

k 

c 

d 

SSFI' 

fT 

Volume 

k'* Y, 

1.8485 

0.6734 

0.6911 

1.8475 

0.0088 

4780.7 

0.9725 

Mean Diameter 

k'* Y, 

1.7435 

0.1157 

0.06859 

-0.6656 

-0.03297 

27.3 

0.9851 

Basal Area 

k'* Y 1 

1.6681 

0.321 

0.053 

0.4103 

-0.01285 

153.9 

0.9702 

Number of Stems 

k'*ln(y,) 

12.26 

0.3139 

0.00505 

2.0566 

0.00975 

27307 

0.9994 
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Test of Models 

Precision indexes 

The relative difference between the prediction value and 
the actual value of stand factors (RE or RD), relative dif¬ 
ference percent ( RDP) and its absolute (IflDPI), standard 
deviation ( SdRDP) and standard error ( SdeRDP) are used 
as precision indexes. The formulae are: 


Table 2.Testing results with dependent data set 


Stand factors 

HE 

RDP 

IRDPI 

SdRDP 

SdeRDP 

Number of 

stems 

0.8142 

0.8205 

0.8207 

0.0333 

0.0972 

Mean 

diameter 

0.1367 

0.1776 

4.583 

0.0932 

0.696 

Basal area 

0.0880 

0.6314 

13.89 

0.1865 

2.7886 

Volume 

0.1516 

2.8266 

12.0875 

0.1534 

1.9806 


RE 


I >7-1 r, 

i _i_ 

Z>7 


x 100 


Testing results with the independent data set 
The results of testing models with independent data set 
are listed in Table 3. 


£ av-K.) 

RDP =-i-z-x 100 ( 14 ) 

n 

f lOVzM 

ZmU y < 

I RDP\ = -i -4-xl00 ( 15 ) 

n 

The relative difference (RE) between the predicted mean 
and the actual mean of stand factors indicates the preci¬ 
sion of the predictor of totals. Percent of relative difference 
indicates the average relative difference between the pre¬ 
dictor and real value for a single stand. Moreover, it reflects 
the difference between the actual value and prediction of 
the total. Thus the RE should be used as a precision index 
of model fit when the model will be used for macro deci¬ 
sion-making (Forest planning). The absolute value of rela¬ 
tive difference percent reflects the average difference of 
actual value and prediction for one stand in the total, so the 
RDP should be used as the precision index when the 
model will be used for micro decision-making (For example, 
forest management). 

Testing results with the dependent data set 
The results of testing models with the dependent data 
set are listed in Table 2. Table 2 shows that the estimates 
of total means of all stand factors have high precision (er¬ 
ror less than 1%). By comparison, the estimates for single 
stand have lower precision. For example, when the total 
volume is estimated, its precision is about 99.85% 
(100-0.1516). But when one stand volume is estimated, its 
precision is about 87.9% (100-12.0875). 


Table 3. Testing results with independent data set 


Stand factors 

RE 

RDP 

IRDPI 

SdRDP 

SDeRDP 

Number of 

stem 

-0.6587 

-0.6973 

0.7007 

0.0702 

0.1515 

Mean 

diameter 

-0.5043 

-1.0764 

5.1633 

0.2206 

1.3558 

Basal area 

-1.312 

-2.7789 

9.597 

0.283 

2.4386 

Volume 

-4.514 

-3.6682 

11.045 

0.3014 

2.9329 


Table 3 shows that the models have a little lower preci¬ 
sion when they are used for an independent data set, but 
the precision are still high. Their RES are less than 5%, 
and \RDP\ are not significantly different from that of the 
dependent data set. 

Coordinating the relationship among the number of stems, 
mean diameter and basal area 

The relationship among the number of stems, diameter 
and basal area is known as following: 

1A 3.1416 

Basal Area=-x number of stem x 

4 

(Mean diameter) 2 (16) 

It is obvious that equation (16) does not exist when the 
basal area, number of stems, and mean diameter are 
separately calculated with models. The precision of the 
basal area calculated by the two methods are given in Ta¬ 
ble 4. It suggests that the basal area can be calculated by 
either of the methods. 


Table 4. Precision comparison of basal area 


Type of data set 

Estimation method of basal area 

RE 

RDP 

\RDP\ 

SdRDP 

SdeRDP 

Independent 

Method A 

-1.3124 

-2.7789 

9.5970 

0.2830 

2.4386 

Method B 

2.9514 

-1.1380 

11.7335 

0.3390 

3.1484 

Dependent 

Method A 

-0.08796 

-0.6314 

13.8900 

0.1865 

2.7886 

Method B 

-.5107 

-2.6603 

10.8900 

0.1428 

1.7177 


Note: Method A: Calculated with mean diameter and number of stems; Method B: Calculated with model directly. 
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Use of models 

The test results indicate that the dynamic stand model 
has high precision and can be used for predicting total or 
single stand growth. Since the length of the growth period 
is 3 years in the original data, it has the most precise pre¬ 
dictors when A2=A1 +3 in practice. 

The state and the dynamic changes in the period of pre¬ 
diction of specific stand are taken into account in dynamic 
stand growth models. Since highly accurate predictors can 
be obtained with the models, the models can be used for 
stand management strategies. 

The model can be used for forest resource file man¬ 
agement. Their highly precise predictors can be used for 
updating the data of the compartment in resource file man¬ 
agement. 

Dynamic stand models can be used for predicting the 
stand diameter distribution and timber ratio. Substituting 
the prediction of the models for the parameter prediction 
models or parameter recovery models to calculate the pa¬ 
rameters of stand diameter distribution and estimate the 
number of stems for every diameter class further. The tim¬ 
ber ratio can then be estimated with the assortment model. 

Conclusion 

The dynamic stand models for larch plantations were 
established with continuous forest inventory data. The re¬ 
sults of tests indicated that the models had high precision. 
They also showed that the continuous forest inventory data 
can be used not only for macroscopic decision-making, but 
also for stand management strategies. 

The stand site quality had to be evaluated with average 
stand height because of the absence of dominant height 
data. But the results of testing the models show that site 


class index can evaluate site quality. However, even better 
results can be expected with suitable site index tables. 

Accurate predictions for specific stands produced with 
the models can be used for stand management strategies 
and forest resource management, especially for renewing 
compartment data in forest resource files. 
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